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Abstract. 

We measured high resolution Cu L3 edge resonant inelastic x-ray scattering 
(RIXS) of the undoped cuprates La2Cu04, Sr2Cu02Cl2, CaCu02 and 
NdBa2Cu30e. The dominant spectral features were assigned to dd excitations 
and we extensively studied their polarization and scattering geometry dependence. 
In a pure ionic picture, we calculated the theoretical cross sections for those 
excitations and used them to fit the experimental data with excellent agreement. 
By doing so, we were able to determine the energy and symmetry of Cu-3d states 
for the four systems with unprecedented accuracy and confidence. The values of 
the effective parameters could be obtained for the single ion crystal field model but 
not for a simple two-dimensional cluster model. The firm experimental assessment 
of dd excitation energies carries important consequences for the physics of high T c 
superconductors. On one hand, having found that the minimum energy of orbital 
excitation is always > 1.4 eV, i.e., well above the mid-infrared spectral range, 
leaves to magnetic excitations (up to 300 meV) a major role in Cooper pairing in 
cuprates. On the other hand, it has become possible to study quantitatively the 
effective influence of dd excitations on the superconducting gap in cuprates. 
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1. Introduction 

Despite the enormous effort devoted to the study of high T c superconductivity in 
cuprates, a general consensus on the underlying mechanisms is still lacking. It is 
commonly agreed that low energy elementary excitations should play a crucial role 
in the formation of conduction electron Cooper pairs, so that considerable efforts 
are being devoted to make a link between phonon and magnon spectra and high 
T c in cuprates pQ. Lattice modes (phonons), which are at the basis of conventional 
superconductivity as explained by BCS theory, lie below 90 meV, an energy seemingly 
incompatible with critical temperatures peaking as high as 130 K. And the possible 
role of spin excitations (magnons) is still debated: the superexchange coupling in 
layered cuprates is exceptionally large (J ~ 130 meV) and gives rise to magnetic 
excitations up to almost 300 meV. Although in doped superconducting materials the 
long range magnetic order is lost, short range magnetic correlations persist [2], and 
magnetic excitations survive across the whole reciprocal space as recently shown by 
Braicovich et al [3J. The possible role of magnetic excitations onto Cooper pairing 
has been the object of several works in the past [4] and is still very actively studied 
at present [5]. 

In this context orbital (dd) excitations, which correspond to a change in the 
symmetry of the occupied Cu-3d orbitals, have been attracting much less attention, 
because their energy is generally much higher (typically > 1.5 eV). However, 
starting from an advanced use of Eliashberg equations [B], Little et al [7] argued 
that dd excitations are also possibly implicated in the mechanisms of high T c 
superconductivity. Indeed a relation between the Cu to apical-oxygen distance and T c 
had been found already several years ago by Ohta and coworkers [5], who proposed 
a rationale based on Madelung potentials. They summarized their results in the so- 
called "Maekawa plots", reporting the dependence of T c on the energies needed to 
transfer the Cu-3d a .2_ y 2 hole to a Cu-3<i3 2 2_ r 2 orbital or to the 0-2p states. More 
recently Sakakibara et al |9j have refined the theoretical analysis by a two-orbital 
model applied to the model structure of La2-a;Sr x Cu04 and HgBa2Cu04+5. In their 
calculations the key parameter is the energy difference between the Cu-3d x 2 _ y i and 
Cu-3d 3z 2_ I ,2 states (AE eg ): d-w&ve superconductivity is favored by a larger energy 
splitting because of a reduced contribution of the d 3z 2_ r 2 state to Fermi surface. In 
the mentioned works the authors based their discussion on theoretically calculated dd 
excitation energies rather than on experimental ones. However, there is no consensus 
on the actual values because theoretical and experimental estimates vary from author 
to author. It is thus timely to provide to the theory of superconductivity a firmer 
experimental basis. Here we provide all the dd excitation energies in several undoped 
layered cuprates measured by resonant inelastic x-ray scattering (RIXS) at the Cu L3 
edge. 

The Cu ion in 2D layered cuprates is nominally divalent (Cu 2+ ), corresponding 
to a 3d 9 electronic configuration. It is generally accepted that in the ground state 
the 3d hole on Cu is mainly found in a d x i_ y i orbital [10] . dd excitations correspond 
to a change in the symmetry of occupied 3d states. In particular, for Cu 2+ , the 
hole is excited from the d x 2_ y 2 symmetry to the d3 2 2_ r 2, d xy , d xz or d yz orbitals. 
dd excitations are an important ingredient for the understanding of the electronic 
structure of cuprates beside their possible direct implication in the mechanism of high 
T c superconductivity itself. In analogy we can look at NiO, whose dd excitations 
have been studied by optical spectroscopies [11] [12] , electron energy loss spectroscopy 
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(EELS) [HE] and, recently, by RIXS [J5] [TBI El- However, although extensively 
investigated, dd excitations in cuprates are still object of debate because of the lack 
of conclusive experimental results. 

Most of the experimental basis comes from optical spectroscopy. In the optical 
absorption of I^CuCU and Sr2Cu02Cl2, Perkins et al [18] found a sharp feature at 
0.41 (0.36) eV, which they assigned to transitions to the dy, z 2_ r 2 orbital. The same 
authors identified the transition to the d xy orbital to occur at 1.50 eV in Sr2Cu02Cl2. 
Elcctrorefiectance measurements by Falck et al [19] on La2CuC>4 revealed the 
transitions to the d xy and d xz / yz orbitals to be at 1.40 and 1.60 eV, respectively. This 
scenario was compatible with crystal field calculations accompanying Cu L% RIXS 
measurements by Ghiringhelli et al [20] that located the transitions to the d 3z 2_ r 2, 
d xy and the doubly degenerate d xz / yz orbitals at 0.41 (1.17), 1.38 (1-29) and 1.51 
(1.69) eV for La 2 Cu04 (S^CuC^C^). However, Lorenzana and Sawatzky [3TJ had 
previously proposed that the feature at ~ 0.4 eV in the optical absorption spectra is 
rather due to a phonon assisted bimagnon excitation. And more recent RIXS results 
by Braicovich et al [22] confirmed the latter assignment for La2Cu04 and CaCuO^: 
the mid-infrared feature is due to magnetic excitations, dispersing up to ~ 0.40 eV, 
and strong dd excitations are found around 1.5 - 2.5 eV. This fact is also compatible 
with the optical Raman scattering measurements by Salomon et al |23j , who identified 
the transition to the d xy orbital at 1.70 eV for La2Cu04 and 1.35 eV for S^CuC^C^. 
Moreover Cu edge RIXS measurements in S^CuC^C^ by Kuiper et al [24] had 
already located the transitions to the d 3z 2_ r 2, d xy and d xz / yz orbitals at 1.50, 1.35 
and 1.70 eV, respectively (see also [20] for Cu L 3 RIXS results). Finally direct first 
principle calculations by Middlemiss et al [25] support these findings too. 

By taking advantage of the recent experimental improvements in the field of soft 
x-ray RIXS we address here the problem of determining the dd excitation energies in 
cuprates in a systematic way. We have measured Cu L$ RIXS spectra of La2Cu04 
(LCO), Sr 2 Cu0 2 Cl 2 (SCOC), CaCu0 2 (CCO) and NdBa 2 Cu 3 6 (NdBCO) at several 
sample orientations, dd excitations are shown here to have local character with no 
or very little dispersion vs the in-plane transferred momentum q» . Using single ion 
theoretical cross sections we could determine the dd excitation energies for all samples 
with a high degree of confidence. 

2. Resonant inelastic x-ray scattering 

-^2,3 edge RIXS is emerging as a powerful technique for the study of neutral excitations 
in 3c? transition metal (TM) oxides and cuprates in particular. Tanaka and Kotani |26] 
suggested that dd and charge transfer (CT) excitations in cuprates can be seen with 
this technique. It should be emphasized that, while forbidden in the optical absorption 
spectra by dipole selection rules, dd excitations are allowed in RIXS due to two 
consecutive dipole transitions. In fact, the scattering process taking place in L 3 edge 
RIXS can be seen as follows. Initially a 2p 3 / 2 electron is resonantly transferred in the 
3d shell through the absorption of an x-ray photon: the system is then in a highly 
excited state with a deep core hole. Secondly the system decays via the transition of 
a 3d electron into the 2p 3 / 2 states and the emission of a photon. As the intermediate 
state, characterized by a 2p core hole and an extra 3d electron, is not observed, the 
whole process has to be described at the second order by the Kramers-Hciscnbcrg 
equation. We deal then with an energy loss spectroscopy, i.e., an inelastic scattering 
of x-ray photons that leaves the solid in an excited state. The energy and momentum 
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of the final state are known from the measured variation of energy and momentum of 
the scattered photons. In this work we restrict to the cases where the system is left 
in a final configuration corresponding to a dd or spin-flip excitation. 

It must be noted that local spin-flip excitations are not eigenstates of the 2D 
antifcrromagnetically ordered lattice. Rather spin waves (magnons) are excited by Cu 
L3 RIXS. Those collective excitations are known to disperse in energy as function of 
their momentum and are traditionally mapped by inelastic neutron scattering (INS). 
In the specific case of layered cuprates magnetic excitations are particularly difficult 
to measure with neutrons due to their very high energy and only rather recently 
high quality results were obtained for La2Cu04 [47]. Very recently the same type of 
measurements were made by Braicovich et al using Cu L 3 RIXS [22[ [3[ [28] . Thanks to 
the relatively large momentum carried by soft x rays, the q-dependence of single and 
multiple magnons could be investigated over approximately 2/3 of the first Brillouin 
Zone (BZ). Those results were also supported theoretically [37] ■ It is thus clear that, 
thanks to the advances in the instrumentation [3T] [33] [3U] that led the experimental 
linewidth in the 100 meV range, Cu L3 edge RIXS has become a complementary 
technique to INS for the measurement of spin waves dispersion. 

On the other hand dd excitations have been studied by Cu L3 [20] and M2,3 edge 
(3p-3d) RIXS [24] with partial success. At M the lower photon energy provides a 
potentially better energy resolution as compared to L edge; however, because of the 
limited momentum carried by photons the portion of the BZ which can be probed 
is limited to a region close to the F point. Any possible dispersion is thus hardly 
detectable. Furthermore the insufficient spin-orbit splitting of the 3p states and the 
extreme weakness of the inelastic cross section with respect to the elastic scattering 
make the analysis and interpretation of the experimental data very difficult. It is 
interesting to note that dd excitations have also been seen at Cu K edge RIXS [33] 
but in this case dd are much weaker than CT excitations because they are indirectly 
excited. 

3. Experimental information 

LCO, SCOC and NdBCO measurements were performed at the ADRESS bcamline [29] 
of the Swiss Light Source at the Paul Scherrer Institute (SLS-PSI) in Villigcn (CH) 
using the SAXES [3D] spectrometer. At the Cu L 3 edge (approximately 930 eV) the 
combined energy resolution is AE = 130 meV FWHM with a data point sampling 
of 14.1 meV/data-point. CCO RIXS spectra were recorded at the ID08 beamline 
of the European Synchrotron Radiation Facility (ESRF) in Grenoble (FR) using the 
AXES [31] [32] spectrometer. The combined energy resolution here was AE = 240 
meV FWHM (31.0 meV/data point) at the same incident energy. In both cases, 
the incident photon energy was finely tuned at the L3 peak by inspection of the 
absorption spectrum. The linear polarization of the incident beam could be set either 
perpendicular (vertical polarization, a, hereafter) or parallel (horizontal, tt) to the 
scattering plane. The polarization of the outgoing beam was not detected. Spectra 
were obtained as the sum of 6 to 12 partial spectra of 5 or 10 minutes each. Off 
resonance elastic peaks were periodically measured on a graphite powder to determine 
the zero energy loss of the spectra. The pressure in the measurement vacuum chamber 
was better than 3 x 10 -9 mbar. All the measurements were done at 15 K (except LCO 
that was measured at room temperature). 

In figure [1] the scattering geometry is shown. The beam hits the sample surface 
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Figure 1. Color on-line. Panels a and b: The experimental geometry is shown. 
The incoming beam hits the sample surface (assumed to be parallel to the ab 
plane) at incident angles 8i and <f>i an( i the outgoing beam is collected at the 
angles o and <j> . The scattering angle 28 is fixed, while the incident angle and 
the azimuthal angle can be changed. They define (5, the angle between the sample 
c-axis and the transferred momentum q (red arrow). The projection of q onto the 
sample ab-plane, , is also shown. In the experiment <5 is changed by rotating the 
sample around a vertical axis. In this way the regions of the 2D reciprocal space 
indicated by thick lines can be covered (panel a, to the right). The [1,0] and [1,1] 
directions correspond to <j>i = and 45° respectively. Panel c: the orientation of 
the atomic spin moments is also a parameter in the cross section calculations. 



at incident angles 9i and <pi and the outgoing beam is collected at angles 9 and <p 
in the xyz reference system; in general the scattering angle 29 is determined by 9i, 
<pi, 9 and <p , but in our case <p a = 4>i so that 29 depends only on 9i and 9 . The 
momentum q transferred to the sample due to the scattering process is shown by 
the red arrow. Its projection onto the sample aft-plane (qu) is also shown since the 
meaningful reciprocal space is usually 2D. The scattering angle 29 is fixed at 130° for 
the AXES spectrometer, while it could be set either at 90° or 130° for SAXES, in 
all cases lying in the horizontal plane. It should be noted that, for a fixed 29, q is 
fixed, but q« can be easily changed both in magnitude and direction in the sample 2D 
reciprocal space by rotating the sample itself around an axis perpendicular or parallel 
to the scattering plane. For the sake of simplicity we introduce the angle 5 between 
the transferred momentum q and the sample c-axis, laying in the scattering plane. 
With this notation, g.i = 2 |fe| sin(0) sin(5), where k = k{ (« k a ) is the momentum 

carried by ingoing (outgoing) photons. For example, when 5 = (specular geometry), 
q is parallel to the sample c-axis and <Z||=0 thus allowing us to probe excitations 
in the center of the BZ. On the other hand, if one could go at infinitely grazing 
incidence (emission), then 5 = —9 (+9) and the maximum magnitude of qu would 
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then be reached. The value of q depends on the scattering angle: at hxii = 930 eV 
q|| mox = 0.47 A" 1 for 29 = 90°, and 0.77 A" 1 for 29 = 130°, respectively. All the 
measurements were performed at a fixed scattering angle 28 and changing S by rotating 
the sample around an axis vertical in the laboratory space and perpendicular to the 
scattering plane, in steps of 5° typically. Samples were aligned so to have parallel 
to the [10] direction (fa = 0) or the [11] direction (fa = 45°) in the BZ. 

LCO and CCO 100 nm thick films were grown by pulsed laser deposition on 
(001) SrTiOs single crystals. Excimer laser charged with KrF (A = 248 nm, 25 
ns pulse width) was used at laser fiuence at about 2 J /cm 2 at the target. Growth 
temperatures were 800 °C for LCO and 650 °C for CCO. A partial oxygen pressure 
of about 0.1 mbar was employed to correctly oxidize the films during deposition. In 
the case of LCO the quality of the films was also in-situ controlled by reflection high- 
energy electron diffraction (RHEED) technique. The SCOC sample were laminar-like 
single crystals obtained by slow cooling in air of a melt of Sr2Cu02Cl2 powder. X-ray 
powder diffraction (XRD) analysis was conducted on a Rigaku x-ray diffractomctcr 
with Cu Ka radiation (A = 1.5418 A) and the electron microprobe was used for 
chemical analysis . Ndi^Bai.sCuOe+x (x < 0.1) 100 nm thin film were deposited 
on SrTiOs (100) single crystals by diode high-pressure oxygen sputtering. Undopcd 
NdBCO thin films were obtained by reducing the oxygen content by annealing as 
grown Ndi.2Bai.sCu07_,5 samples in Argon atmosphere (10 mbar) for 24 hours. The 
structural and morphological properties of the sample have been checked by x-ray 
diffraction, also using synchrotron radiation, and by atomic force microscopy |35| . 
The samples are tetragonal and perfectly matched with the STO lattice ( a=6=3.905 
A) while the c-axis is 11.81 A. 

For LCO, CCO, NdBCO and SCOC the fraction of the BZ that could be 
explored along the [10] direction was 93.1%, 94.4%, 95.6% and 97.3%, respectively, for 
20 = 130°. In the case of LCO and SCOC, it decreased at 56.9 (59.4)% with 29 set at 
90°. Moreover in the case of SCOC, measurements along the [11] direction were taken 
up to 68.8% of the nuclear BZ boundary with 29 = 130°, i.e., the whole magnetic BZ 
could be spanned. 

4. Experimental spectra and cross section calculations 

4-1- Overview of experimental spectra 

In the left panel of figure [2] a Cu L3 edge RIXS spectrum with a polarization 
of LCO is shown. The scattering angle 29 was set at 130° and the sample was 
rotated at 6 = +45°, i.e., 20° grazing emission, and fa = (q« parallel to the [10] 
direction in reciprocal space). The very high resolution allows us to clearly recognize 
different excitations depending on the energy scale. The spectrum is dominated by 
dd excitations, which will be extensively discussed below. Both at higher and lower 
energies other features are found with a much (approximately one order of magnitude) 
lower counting rate. Excitations in the mid-infrared region (up to ~500 meV) are 
expanded in the inset and are known to have a magnetic character: they are a 
combination of a dispersing magnon, a continuum given by bimagnons and other 
multi-magnon excitations, phonon peaks and an elastic zero- loss line (22j El [28] ■ At 
higher energies CT excitations are also visible, as a broad distribution, via the Cu 
3d-0 2p hybridization. By looking at the absorption spectrum (dashed line), it can 
also be noted that features above 1 eV energy loss, such as dd excitations, are weakly 
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Figure 2. (Left panel) Example of Cu L3 absorption (dashed) and RIXS (solid) 
spectra of LCO with a polarization. For RIXS the angles were set to 2f5 = 130° , 
8 = +45° (i.e. 20° grazing emission) and (f>i = 0. One can immediately recognize 
charge transfer (CT), dd and magnetic excitations at different energy ranges. 
A closer look of the mid-infrared energy region is given in the inset. (Right 
panel) RIXS spectra for LCO, SCOC, CCO and NdBCO in the same experimental 
geometry as left panel. 

affected by self-absorption because of the reduced absorption coefficient 2 eV below the 
resonance peak. In the right panel of figure [2 dd excitations of LCO are compared to 
those of SCOC, CCO and NdBCO both for a (solid) and 7r (dashed line) polarization. 
The effect of changing incident photon polarization is considerable. The spectral 
shape changes drastically for the four compounds, reflecting differences in the Cu 2+ 
coordination. In order to recognize the symmetry of the RIXS final state we carried out 
systematic measurements with both polarization as a function of 6 (or 8i, accordingly) 
for a fixed 29 and fitted the spectra to obtain an estimate of the dd excitations energies. 
The fitting procedure is presented in the following sub-section and is based on single 
ion model calculations of the RIXS cross sections. 

4-2. Single ion model cross section calculations 

In cuprates, Cu ions are known to be mostly in the Cu 2+ oxidation state [39] 
corresponding to a 3d 9 electronic configuration. In the crudest approximation, wc 
only consider the atomic states of a Cu 2+ ion. RIXS is a second order process and is 
described by the Kramers-Hciscnbcrg (KH) equation. The scattering is here modeled 
in two steps: first one 2p 3 / 2 electron is resonantly promoted into the 3d states by 
absorption of a photon (intermediate state). Because of the large 2p spin-orbit splitting 
(~ 20 eV), interference effects with the 2p 1 / 2 states can be neglected. The only 
available intermediate state is given by a fully occupied 3d shell and one hole in the 
four- fold degenerate 2p 3 / 2 core level. The second step is the radiative deexcitation of 
one 3d electron into the 2p 3 / 2 levels. In a short notation, the process can be written as 
the sequence 2pj /2 3d 9 -> 2p\ /2 3d w ->• 2pj /2 3d 9 *, where the * indicates that the final 
states can be either the ground state or an excited state with the 3d hole occupying a 
state with different orbital symmetry and/or spin. 

Matrix elements entering the KH equation for the calculation of the RIXS cross 
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Figure 3. Single ion Cu 2+ L 3 edge RIXS cross section for 29 = 90° (left 
panel) and 26 = 130° (right panel) and for two different angles (pi = and 
<pi = 45° both for cr and it polarization. The bottom panels show sketches of 
the experimental geometry. The spin is fixed in the calculation along the [110] 
direction (9 S = 90°, (f> s = 45°). All the possible dd excitations including the 
elastic and spin-flip final states are considered. As already explained in figure [l] 
to a given S (bottom axis scale) corresponds an in-plane transferred momentum 
9 II (top axis scale). In the top panels the sums of all dd excitations (excluding 
elastic and spin-flip) cross sections are given. 
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sections are here calculated in the atomic approximation as follows. We consider the 
scattering ion as a hydrogen-like system of a single positively charged particle (hole) 
and we express the angular part of the one-particle wave functions (p and d orbitals) in 
spherical harmonics. The radial part of the integrals is the same for all the considered 
transitions and is thus neglected also in the expressions of the 2p and 3d states. The 
local symmetry is , so the commonly used set of atomic d orbitals are well suited 
to represent the atomic states in terms of symmetry: 

d X 2-y2 = ~J= (F 22 + F 25 ) (1) 
d 3z 2„ ? ,2 = Y 20 (2) 

d xv = - ^= (Y22 - Y 22 ) (3) 
d yz = --j=(Y 21 +Y 2l ) (4) 

d xz = ^=(Y 21 -Y 2l ) (5) 

The latter two are degenerate for symmetry. As in all layered cuprates the four in- 
plane O ions are the Cu nearest neighbors, the unique 3d hole present in the ground 
state has an almost pure x 2 — y 2 character [10] . Thus the dd excitations can be 
described as the transfer of the 3d hole from the d x 2_ y 2 orbital to the d 3z 2_ r 2, d xy , d yz 
or d xz orbitals. Also the case of d x 2_ y 2 final state can be viewed as a dd excitation, 
giving rise to an elastic scattering or to a pure spin-flip excitation depending on the 
final spin state. 

Superimposed to the CF splitting of 3c? orbitals is the additional spin splitting 
of each state due to super-exchange interaction with neighboring in-plane Cu 2+ ions 
(spin orbit of 3d states is neglected) which doubles the number of possible final states. 
In this work we will not consider the low energy part of the spectra because the 
description of the set of possible final states cannot be done in a single ion model. 
In fact, the local spin-flip is not an eigenstate of the 2D Heisenberg AFM lattice: 
spin waves should then be introduced to account for the dispersion of pure magnetic 
excitations, as already demonstrated in references [2"T1 13"]. 

Although we do not look at pure magnetic excitations, in order to correctly 
calculate the RIXS cross sections the local atomic spin orientation has to be explicitly 
taken into account. As the 3d spin-orbit interaction is weak in Cu [34] spin is a good 
quantum number for both the ground and the final states. We use here the Pauli 
matrices & x ,a y ,a z . The arbitrary spin orientation of the initial hole is defined by the 
angles 6 S and cj> s and the eigenvectors of 

<j(0 s , 4> s ) = (<7 X cos 4> s + a y sin <j> s ) sin 6 S + a z cos 9 S (6) 

give the weights for spin up and spin down components along z for the generic spin 
direction. For instance, if we assume the ground state spin direction to be "down" in 
the hole representation, the ground state wave function and the- spin-flip wave function 
are written as 



d^ 2 2 



d\ 



1 

V2 
1 

V2 



U- (Y 2 \ + 4) + D_ (Y 2 \ + 4) ] (7) 
U+ (Y 2 \ + 4) + D+ (4 + 4)1 (8) 
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where {/_({/+) and D_(D + ) are the components of the eigenvector corresponding 
to the negative (positive) eigenvalue of cf(6 s ,<j) s ). The same holds for the other 
3d orbitals which, together with the ground state itself and the spin-flip state, 
represent the 10 different possible final states available for the transition, i.e 

^t 2 -y 2 ' ^x 2 -y 2 ' ^3z 2 -r 2 ' ^L 2 -r 2 ' ^xy: d\ y , d^ z , d\ z , d^ z and S yz . 

On the contrary, 2p states are strongly spin orbit coupled and spin is not a good 
quantum number; the energetically degenerate 2j> 3 / 2 orbitals are written in terms of 
spherical harmonics as follow: 

Ph-I= Y A ( 9 ) 
P s x = JIy* + J\y& (10) 



= *n (12) 



In the approximation that the lifetimes of the 2p 3 / 2 states are the same, the RIXS 



transition from the ground state d 



to a given final state, say d 



4(t) 



T ±(t) 

7 3z 2 - 



32 2_ r 2, is given by 
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E 



(d^_ r2 \T^\p h J(ps !m \T q \di 2 



(13) 



runs over the intermediate states. 



T„ 



y/4.it/{2q + l)Yig(0, 4>) is the expression for the electric dipole operator in spherical 
harmonics with q — —1,0, 1 for left, linear (z) and right polarized light, respectively. 
The sum over q' is required since we do not measure the polarization of the outgo- 
ing photons in our experiments. Matrix elements are thus simple integrals of three 
spherical harmonics which can be readily calculated [40] . 

Using this method we can calculate the scattering cross section for all the dd 
excitations of cuprates for any incident and scattered directions and polarization of 
the photons and any atomic spin direction. For example, if the atomic spin is oriented 
along the [110] direction (8 S = 90° and 4> s = 45°) and 29 = 90°, the cross sections for 
the elastic and pure spin-flip transitions are given by the simple formulas 
I 

x —y 



cx 4 



(14) 



o~ x 2_ y 2 :Cr oc sin 2 6 l (15) 

in the case of incident a polarization and 

ff^K sin 2 20, (16) 

o-\i_ y 2 >7r oc cos 2 6i (17) 

in the case of incident tt polarization and the scattering plane perpendicular to the 
sample surface. 

In figure [H the calculated cross sections are shown in the two particular cases 
<pi = and <f>i = 45° as a function of 9i or alternatively of the in-plane transferred 
momentum qu for 29 = 90° (left panel) and 29 = 130° (right panel). We underline 
that, contrarily to what previously stated in the literature [301 HI], pure spin-flip 
transitions are allowed as long as the spin is not parallel to the [001] direction [27] 
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(see the curves labeled a f and ir f for the x 2 — y 2 final state in figure [3]). For layered 
cuprates, this selection rule has the important consequence that single magnons can 
contribute to the Cu L 3 edge RIXS spectrum [37J since here spins are known to 
lie always in the a6-plane [HI 23J [35] ■ We highlight the fact that the orientation 
of the spin within the ab plane has no influence on the spin-flip cross section. This 
makes it simpler to measure the magnon dispersion in samples with multiple magnetic 
domains. However, some of the dd cross sections do depend on the <j) s value (having 
fixed 6 a = 90°). By this dependence one could think of using the dd excitation 
spectrum to determine the in-plane orientation of the spin in layered cuprates. 

Calculations are then used to determine the energy and symmetry of Cu-3<i states 
by fitting the experimental data. It is important to note that cross sections are here 
calculated on the basis of symmetry properties of the angular part of the ground 
state and final state wave functions. By including the Cu Zd-0 2p hybridization 
the symmetry of the problem would not change [44 [ |4"5] and the dependence on the 
scattering geometry and photon polarization would remain the one calculated here, as 
it depends solely on symmetry properties [45) . 

Once the cross sections for all the possible final states are calculated, one has to 
construct the simulated spectra as a function of the energy loss (Huj — hui) as 



where / runs over x 2 — y 2 , 3z 2 —r 2 , xy, xz and yz. We assume here that the Lorentzian 
lifetime broadening is equal for states with the same orbital symmetry while the 
super-exchange coupling (Jf) is orbital dependent. From overlap considerations |46j 
one finds that J 32 2_ r 2 = J x 2_ y 2/6 = J/6 and J xy = J xz = J yz = 0. J/ gives the 
energy separation between a peak and its spin-split satellite: according to the 2D 
Ising model, in the ground state the energetic cost to flip one spin is 2 J, while in the 
excited states the exchange constant is reduced because of the smaller overlap with 
nearest neighbor orbitals. In the present work we have fixed J = 130 meV for all 
samples for simplicity, although it is known that its value is different (±20%) from 
sample to sample. Indeed this difference would little affect our calculated spectra, due 
to the strong reduction or cancelation of Jf with respect to J as explained above. 

Finally Ef is the energy of the final state with a given symmetry (/) and spin 
down (J,) . Geometry and polarization dependencies that enter the simulations through 
o~f are obviously fixed by the experiment. This imposes a severe constraint to the 
number of free fitting parameters which is thus limited to the energy positions of the 
3d states (namely E 3z 2_ r 2, E xy and E xz / yz , while E x 2_ y 2 = by definition) and to 
their Lorentzian lifetime broadening (r 3z 2_,,2, T xy and Y xz / yz ). The energies (Ef) of 
the final states should be considered as effective energies, thus taking into account 
both the ionic and the covalent part of the bond. Figure 0] shows an example of the 
fitting procedure for LCO (6 = 0, 26 = 90°, a polarization). Thick blue lines are <5-likc 
functions whose heights are proportional to the calculated cross sections ((Jf); they 
are convoluted with Lorentzian functions to take into account the lifetime broadening 
of the final states (r^) to obtain the dashed curves. Their sum gives the red dashed 
line which is eventually convoluted with a Gaussian (whose FWHM matched to the 
effective experimental resolution) to obtain the simulated Cu L 3 edge RIXS spectrum 
(red line). 



I(Huj — hijji) 



(18) 
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Figure 4. Cu L-$ edge RIXS spectrum (open circles) of LCO taken at 5 = 0, 
29 = 90° and a polarization. The theoretical spectrum is built by convoluting 
delta functions proportional to the calculated cross sections a f (thick blue) with 
Lorcntians to include finite lifetime broadening (dashed blue). Their sum gives 
the hypothetical spectrum measured with infinite resolution (dashed red). The 
spectrum to be compared to the experimental results is obtained after convolution 
with a Gaussian curve (solid red). The energy of the various dd excitations (Ef) 
and the Lorentian widths (T f) are used as adjustable parameters in the fitting 
procedure. 



5. Determination of dd excitation energies and discussion 

5.1. Fitting results 

In this section, the experimental trends are shown together with the results of the 
simulations. In figure[5l measured and calculated RIXS spectra of LCO are compared, 
both for 29 = 90° and 29 = 130° with 0j = 0. From the top spectrum to the bottom 
one 5 is changed in steps of 5° from 10° grazing emission to 10° grazing incidence. As 
the incidence angle is changed, q^ varies accordingly. However, it is evident from the 
experimental results that the dispersion of dd excitations, if any, is very small both 
for 29 = 90° and 29 — 130°. This finding seems to be general and holds true also in 
the cases of SCOC (figures M and [7]), CCO (figurcHJ) and NdBCO (figurci]). Within 
the present experimental accuracy the dd excitations do not disperse in any of the 
undoped layered cuprates measured by us. 

The fact that dd excitations are fixed in energy when changing S allowed us to 
fit all the spectra for a given sample using fixed values of Ef and Tf. Using the 
redundance of the experimental basis the values of the free parameters Ef and Tf 
could thus be optimized. The high degree of confidence reached is made evident by 
the excellent agreement of the fitted spectra to the experimental ones shown in the 
figures. For clarity we remind here the input parameters of our fitting. The RIXS cross 
sections 07 are calculated for each final state within the ionic model; the superexchangc 
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Figure 5. Measured (black line) and calculated (red line) RIXS spectra of 
LCO for various <5, ranging from 10° grazing incidence to 10° grazing emission 
with respect to the sample surface in steps of 5°. Dispersion along [10] direction: 
fa = 0; 29 = 90° (left panel) and 29 = 130° (right panel). 




Figure 6. Measured (black line) and calculated (red line) RIXS spectra of 
SCOC for various <5, ranging from 10° grazing incidence to 10° grazing emission 
with respect to the sample surface in steps of 5°. Dispersion along [10] direction: 
fa = 0; 29 = 90° (left panel) and 29 = 130° (right panel). 



interaction is fixed to J = 130 meV for all samples; the Gaussian broadening is set to 
130 meV FWHM; the energy Ef and Lorentzian broadening T f are the free parameters 
for optimizing the fitting and are thus the result of the comparison of calculated and 
measured spectra. 

The results for the four samples are summarized in Table [TJ For LCO and SCOC 
we used data measured both at 29 = 90°, 130° for fa = 0. For SCOC we employed data 
with 29 = 130° and fa = 45°. For CCO and NdBCO the used data were at 29 = 130° 
and fa = only. Although the spectra of CCO shown in figure [8] were measured with 
lower resolving power (AE = 240 meV), the fitting results are perfectly compatible 
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Figure 7. Measured (black line) and calculated (red line) RIXS spectra of SCOC 
for various <5, ranging from 10° grazing incidence to 10° grazing emission with 
respect to the sample surface in steps of 5°. Dispersion along [11] direction: 
4>i = 45°; 20 = 130°. 

with the (few) spectra taken at higher resolution (AE = 130 meV, figure[2]). However, 
here the better resolution makes evident an additional feature at about 2.4 eV, 
which could not be assigned within our model. Possibly the extra peak could be 
related to oxygen vacancies affecting the local environment of Cu ions. Independent 
results on magnetic excitations of the same sample |49j seem to further support this 
interpretation. 

Finally we notice that the fitting works well for NdBCO too, in spite of the 
presence in this sample of two different Cu sites. Namely in this YBCO-like system Cu 
ions can be found also outside the CuC>2 planes, in the so-called CuO chains (although 
in this undoped compound chains are actually broken). As in the simulations we 
assume that only one specie is contributing to the RIXS spectrum, we attribute this 
result to the extreme selectivity provided by the resonance in the excitation step. 

5.2. Discussion 

The energies of the 3c? states obtained through our RIXS measurements are evidently 
related to the local coordination of Cu ions (symmetry and atomic distances). In 
fact, the 4 materials share the four-fold planar coordination that gives the famous 
CUO4 plaquettes present in all cuprate superconductors. However, the in-plane Cu-0 
distances vary from a minimum in LCO to a maximum in SCOC: the lattice parameters 
are a = 3.80 A for LCO [36], 3.85 A for CCO [37], 3.90 A for NdBCO [38] and 3.97 A 
for SCOC [36]. The out-of-plane coordination is even more diverse. CCO has no 
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Figure 8. Measured (black line) and calculated (red line) RIXS spectra of CCO 
for various <5, ranging from 10° grazing incidence to 10° grazing emission with 
respect to the sample surface in steps of 5°. Dispersion along [10] direction: 
<k = 0; 29 = 130°. 



Table 1. Parameters used in the calculations to fit the experimental data. Also 
listed are the values of the effective CF parameters and those for a pure covalcnt 
picture within a 2D cluster. In all except the LCO case, covalent parameters could 
not be defined (n.d.) as explained in the text. For all samples the superexchangc 
J was fixed to 130 meV. 





La2Cu04 


Sr 2 Cu0 2 Cl 2 


CaCu0 2 


Sro. 5 Cao.5Cu0 2 


NdBa 2 Cu 3 7 


E-iz 2 -r 2 (r3 z 2_ r 2) [eV] 


1.70 (.14) 


1.97 (.10) 


2.65 (.12) 


2.66 


1.98 (0.18) 


E xy (T xy ) [eV] 


1.80 (.10) 


1.50 (.08) 


1.64 (.09) 


1.56 


1.52 (0.10) 


E xz / yz i^xz/yz) [eV] 


2.12 (.14) 


1.84 (.10) 


1.95 (.12) 


1.93 


1.75 (0.12) 




1.70 


1.97 


2.65 


2.66 


1.98 


AE t2g 


0.32 


0.33 


0.31 


0.36 


0.23 


lODq [eV] 


1.80 


1.50 


1.64 


1.56 


1.52 


D s [eV] 


0.29 


0.33 


0.42 


0.43 


0.32 


D t [eV] 


0.11 


0.13 


0.19 


0.19 


0.14 


A pd [eV] 


2.20 


n.d. 


n.d. 


n.d. 


n.d. 


T pd [eV] 


3.20 


n.d. 


n.d. 


n.d. 


n.d. 


T pp [cV] 


0.81 


n.d. 


n.d. 


n.d. 


n.d. 
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Figure 9. Measured (black line) and calculated (red line) RIXS spectra of 
NdBCO for various 6, ranging from 10° grazing incidence to 10° grazing emission 
with respect to the sample surface in steps of 5°. Dispersion along [10] direction: 
6i = 0:29 = 130°. 



apical ligands and, for that, it is often called infinite layer compound. Both LCO and 
SCOC have two apical ligands, symmetric with respect to the basal plane, at larger 
distance than the in-plane O (tetragonally distorted octahedral symmetry); in LCO 
apical O 2- ions are at 2.43 A from Cu 2+ , in SCOC Cl~ ions are at 2.86 A. Finally 
NdBCO has a double layer YBCO-like structure, i.e., only one apical oxygen at 2.27 A 
from the C11O2 plane. 

The relation of the dd excitation energies to the local structural properties of 
our samples is evident and relatively simple. Before comparing the results to any 
theoretical model we highlight these relations in a purely phenomenological way. The 
experimentally determined energies for each final state symmetry are listed in Table Q] 
and are referred to the x 2 — y 2 ground state. A graphical presentation is made in 
figure [POl We use also the results of Sro.5Cao.sCu02 (SCCO, same structure as CCO 
but with larger in-plane lattice parameter a = 3.90 A) are added, although the spectra 
are not shown here |49j . The xy excitation energy increases when decreasing a. When 
displayed in a logarithmic scale a power relation can be highlighted: a best fit to 
our data gives E xy oc a~ i2 . The splitting of the t2 9 states (A t2g = E xz / yz — E xy ) 
is rather small and almost independent of the out-of-plane lattice parameter I and 
of the presence of apical ligands. On the contrary, the splitting of the e g states 
(A eB = E 3z 2_ r 2 — E x 2_ y 2) varies considerably from sample to sample but a univocal 
trend vs the Cu-ligand distances (a/2 and I) cannot be found due to the other 
parameters at play (ligand element and valence, symmetry with respect to the basal 
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Figure 10. (Left panels) The effective CF parameter lODq as a function of the 
lattice constant: linear (lower) and logarithmic (upper panel) scale. (Right panel) 
The splitting of the t2 g and e g states as a function of the out-of-plane nearest 
neighbor distance I. 

plane). However, if we do not consider NdBCO, which has a pyramidal coordination, 
a qualitative trend emerges: the e g splitting incrases when the apical ligand is farther 
from the Cu ion. And the energy of the 3z 2 — r 2 state differs as much as almost 1 eV 
in LCO and CCO. This demonstrates that the local coordination can hugely impact 
on the e g splitting. 

The widely known crystal field (CF) model can be used to understand the trends 
in the experimental results. We take a single site, purely ionic picture, where the 
energies of the 3d orbitals is determined by the symmetry and strength of the non- 
central Coulomb field produced by the charged ligands nearby the central Cu 2+ ion. 
Beside an addictive constant, the energy levels of d states in the tetragonally distorted 
octahedral symmetry are usually written as [50] 



E%F_ y2 =6Dq + 2D s -D t (19) 

££JL r 2 = 6Dq - 2 As - 6 A (20) 

E™ = - ADq + 2D S - D t (21) 

E°J, yz = -4Dq-D s +m (22) 



where Dq (or lODq), D s and D t are the CF parameters depending on the geometrical 
arrangement of point charges around the Cu 2+ ion. In particular, IQDq gives the 
energy splitting between the e g and ti g orbitals, while 4D S + 5D t and 3D S — 5D t 
give the splitting of e g (AE Sg ) and t 2g (AE t2g ) orbitals, respectively (see figure fTTj) . 
It can be shown [50] that within the CF model lQDq scales with a power law of 
the in-plane lattice parameter a, i.e., lODq ~ a~'\ with n = 5. The value of the 
effective CF parameters lODq, D s and D t were determined by comparison with the 
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Figure 11. Local crystal field and energy level diagram for d orbitals in spherical, 
cubic {Oh) an d tetragonal (D^) symmetry. 



experimental results and are given in Table [TJ The single ion CF model neglects 
completely the Cu-ligand orbital overlaps and some important physical properties are 
totally missing, such as the super-exchange interaction. However, it can capture the 
power law dependence of E xy on a, although the exponent n is overestimated (n = 5 
instead of 4.2). On the other hand, although it is always possible to find a combination 
of lODq, D s and D t compatible with the experimental results, it might happen that 
their values have little physical meaning. In particular lODq, which is widely used in 
various calculations from single ion to cluster or impurity models, has to be regarded 
here as an effective parameter, i.e., the value obtained in CF cannot be applied directly 
to other models. 

On the other hand, one can tentatively try to include the effect of covalency 
with the help of a 2D-cluster calculations: Eskes et al [H] proposed a cluster-model 
calculation of the electronic structure of CuO in which they treat the d—d Coulomb and 
exchange interaction within full atomic multiplet theory and use symmetry-dependent 
Cu 3d-0 2p hybridization to describe photoelectron spectroscopic data. In the model, 
they consider a (CUO4) 6- cluster (plaquette), but neglect apical oxygens. In our case, 
this model seems is only to be partially justified: in fact, we can expect that the 
2D-cluster, which docs not include the apical ligangs, cannot fully account for the 
dependence of A e and A t2 on the out-of-plane Cu-ligand distance I. Within this 
covalent picture, the energy levels of d states are written as 
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(26) 



where A p d is the charge-transfer energy, T p d is the ground state Cu 3<i-0 2p 
hybridization energy and T pp is the nearest-neighbor O 2p-0 2p hybridization energy. 
For LCO, the values of the parameters were obtained by taking the value of A pc i from 
independent experiments (2.2 eV [33]) and by fulfilling two of the three equations 
for T p d and T pp . The values which better reproduce the experimentally determined 
energies of Cu-3d states in LCO are T p( j=3.20 eV and T pp =0.81 cV. However, this 
2D-cluster cannot describe the energy sequence found for SCOC, CCO and NdBCO. 
In fact, in those compounds E 3z 2_ r 2 < E xy whereas £§°I_ r 2 < E^,° v for any value of 
A p d, T p d and T pp , as it can be seen from the equations above. It is anyhow interesting 
to notice that this model predicts a power law for the xy to x 2 — y 2 energy splitting 
with exponent n = 3.5, dictated by the dependence of the overlap integrals on the 
Cu-0 distance [46]. This value is rather close to the experimental one. 

6. Conclusions 

Using Cu L 3 RIXS we have unequivocally determined the energy of dd excitations 
(directly related to the energy of Cu-3d states) in several layered cuprates, parent 
compounds of high T c superconductors. These values were largely unknown until 
now, because other experimental techniques could provide only partial or ambiguous 
results. For that reason the possible role of dd excitations in superconductivity of 
cuprates has been often controversial. Three main outcomes can be extracted from 
the comparison of the quantitative results summarized in Table [1] 

i) The dd excitation average energy is about 1.9 eV and, more importantly, the 
minimum lies above 1.4 eV. This finding resolves speculations made in the literature 
about a possible role of 3z 2 — r 2 states lying as close as 0.5 eV to the x 2 — y 2 
ground state. The absence of dd excitations in the mid-infrared spectral region greatly 
supports the hypothesis that magnetic excitations (up to 250-300 mcV) play a major 
role in Cooper pairing in cuprate superconductors, as one can obtain from the solution 
of Eliashberg equations [5]. However, dd excitations can still have an importance for 
the full description of high T c superconductivity [7]. In particular the fact that E 3z 2_ r 2 
greatly varies from sample to sample, is potentially very important in determining T c 



ii) Apical ligands do have an influence on dd excitations, but a quantitative trend 
cannot be predicted by simple models such as crystal field or 2D-clusters. As a 
consequence more sophisticated calculations |52j are needed in order to reproduce the 
energy position of 3d states in cuprates case by determined experimentally. 

iii) A simple relation of the xy state energy to the in-plane lattice parameter exists. 
It is a power law (E xy oc a~ n ) with exponent n ~ 4.2, i.e., in between the prediction 
of crystal field model (purely electrostatic, n = 5) and that of a "covalent" model 
following Harrison's overlap integrals (n = 3.5). The dependence of overlap integrals 
on the interatomic distances is of great importance in cuprates as it drives the super- 
exchange interaction too, which provides the strong antifcrromagnetic background in 
high T c superconductors (nevertheless one should not forget that overlap integrals are 
also strongly affected by Cu-O-Cu bond angles, a variable not considered in this work). 
Again, although indirectly, a better knowledge of dd excitations can be of great help 
in a full description of superconductivity in cuprates. 
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In conclusion wc have made a solid assessment of the dd excitations in 
cuprates. These results can serve as experimental basis for advanced calculations 
of these same energies. Further measurements, both on parent compounds and on 
superconductors, could look for correlations between the dd excitation spectrum and 
the superconducting properties of the cuprates. 
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